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Nonvanishing Local Scalar Invariants even in VSI Spacetimes 
with all Polynomial Curvature Scalar Invariants Vanishing 
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VSI ('vanishing scalar invariant') spacetimes have zero values for all total scalar contractions of 
all polynomials in the Riemann tensor and its covariant derivatives. However, there are other ways 
of concocting local scalar invariants (nonpolynomial) from the Riemann tensor that need not vanish 
even in VSI spacetimes, such as Cartan invariants. Simple examples are given that reduce to the 
squared amplitude for a linearized monochromatic plane gravitational wave. These nonpolynomial 
local scalar invariants are also evaluated for non-VSI spacetimes such as Schwarzschild and Kerr 
and are estimated near the surface of the earth. Similar invariants are defined for null fluids and for 
electromagnetic fields. 

PACS numbers: PACS 04.20.-q, 04.30.-w, 04.70.Bw, 04.20.Jb 



I. INTRODUCTION 

VSI or 'vanishing scalar invariant' spacetimes have re- 
cently been discussed [J 0, i, i, i, i, 0, 0, H and de- 
fined Q as those in which "all of the scalar invariants 
constructed from the Riemann tensor and its covariant 
derivatives are zero." This definition has the implicit 
assumption that 'constructed' means by "contraction of 
a polynomial in the Riemann tensor and its covariant 
derivatives" [1J. 

However, there are other ways of forming local scalar 
invariants, invariants under all passive diffeomorphisms, 
" ofM. For example 0, [H G3, 
l2l| , there are the Cartan invari- 
ants, which are components of the Riemann tensor and its 
covariant derivatives ('curvature tensors') in a frame suf- 
ficiently determined by such curvature tensors that these 
components are invariant under any remaining freedom 
in the frame. (One might say that these other ways are 
'concoctions,' in contradistinction to the limitations of 
'constructions' as implicitly defined above.) 

Here I shall give simple examples from a null vector 
direction determined by curvature tensors (which can be 
one of the null basis vector directions of a null frame 
given by the Cartan method, so my examples can also be 
obtained by that general method). The first examples I 
shall give are local numerical invariants A/i and A/2 that 
are rather simply defined for generic spacetimes (though 
often are either zero or infinite). These are generically 
nonzero for at least certain classes of VSI spacetimes, 
such as plane gravitational waves. These examples are 
concocted by taking the minimum and the maximum, 
over variations of an auxiliary null vector, of the ratio of 
two scalar monomials each formed by the contraction of 



copies of the Riemann tensor or of its covariant deriva- 
tives and eight copies of the auxiliary null vector. For a 
linearized monochromatic plane gravitational wave, both 
of these local scalar invariants are the same and are equal 
to the sum of the squares of the wave amplitudes for the 
two polarizations, h\ + h 2 x . In this case the average value 
of Ai = A/2 may also be interpreted as eight times the 
number of gravitons per wavelength along the direction 
of the wave and per Planck perpendicular area. 

The values of these invariants are also analyzed in the 
Schwarzschild and Kerr spacetimes, even though these 
are not VSI spacetimes and do have nonzero scalar poly- 
nomial invariants, such as the Kretschmann invariant. In 
these non-VSI spacetimes, the nonpolynomial invariants 
M\ and Ni cannot be interpreted as the sum of squared of 
gravitational wave amplitudes or as the number of gravi- 
tons per wavelength and per perpendicular area. 

Other examples are given that can be nonzero for null 
fluids. A third pair is exhibited that are generically 
nonzero for plane electromagnetic waves and give the 
number of photons per wavelength and per perpendic- 
ular area for monochromatic waves. 



II. FIRST TWO INVARIANTS 

For the first pair of invariants, N\ and A/2, define the 
covariant quantities 



i/(k) = R\ lh R\ ld k a k h k c k\ 



(2.1) 



<5(k) = R\ ]b . iCd R' eifwh k a k b k c k d k e kf^k\ (2.2) 
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which for a given event in a given spacetime are functions 
of an auxiliary null vector k with contravariant compo- 
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nents k a . Then define the k-dependent local fraction 



/(k) = 



W] 2 
2<5(k) 



(2.3) 



Because there are the same number of null vector com- 
ponents in the numerator as in the denominator, the ratio 
/(k) is invariant under rescaling k, so it only depends 
upon the direction of k over the unit sphere of direc- 
tions on the null cone. Assuming that the denominator 
is generically nonzero for an arbitrary k, the fact that 
the unit sphere is compact implies that the ratio has a 
minimum, which is defined to be the value of Mi ■ If (but 
not only if) the denominator is nonzero for all nonzero 
k at a given event, then the ratio also has a maximum, 
which is defined to be the value of M%. That is, one may 
define the two local scalar invariants to be 



Mi = min/(k), Mi = max/(k). 

k k 



(2.4) 



If in D dimensions one goes to a frame in which k is 
one of two null basis vectors and the remaining D — 2 
basis vectors are orthonormal spatial vectors orthogonal 
to both of the null basis vectors, then in both Eqs. ()2.1j) 
and (|2.2|) the nonzero terms in the Einstein summation 
over the indices i and j are restricted to these spatial 
vectors, so that both v{k) and S(k) are sums of perfect 
squares. Since each tensor in the product is symmetric in 
i and j, there are generically (D — 2) (D — l)/2 indepen- 
dent squares in each sum. This is larger than the D — 2 
parameters for the direction of k for D > 4, so generi- 
cally in spacetimes with four or more directions, both the 
numerator and the denominator of the ratio /(k) are pos- 
itive for all nonzero k. This allows both Mi and M2 to be 
defined as finite positive local scalar curvature invariants 
for generic spacetimes of at least four dimensions. 

However, in a spacetime that solves the vacuum Ein- 
stein equations, with the Ricci tensor R a b being either 
zero (the case with no cosmological constant) or propor- 
tional to the metric (when the cosmological constant is 
nonzero), the trace of each tensor in the product in either 
j/(k) or <5(k) is zero when i is set equal to j and summed, 
so then there are only \{D — 2)(D — 1) /2] — 1 independent 
squares in each sum, which equals the D — 2 number of 
directional parameters of k for D = 4. Therefore, in a 
four-dimensional vacuum spacetime, generically both the 
numerator and the denominator of the ratio /(k) can be 
zero for discrete directions of k. Indeed, if k is one of 
the principal null directions (PNDs) of the Weyl tensor, 
which always exist in four dimensions, then the numera- 
tor will vanish [22j . Unless the denominator also happens 
to vanish for the same null direction, one gets M\ = 0, 
the generic case for a vacuum four-dimensional space- 
time. On the other hand, there may or may not be null 
directions for which the denominator vanishes in a vac- 
uum four-dimensional spacetime. Both possibilities are 



generic, since the number of equations to be solved, 2, 
matches the number of parameters for null directions, so 
that one can have generic vacuum spacetimes, or regions 
of a spacetime, in which A/2 is infinite, and other generic 
vacuum spacetimes or regions in which M% is finite. 

In a higher-dimensional vacuum spacetime [22I ] , the nu- 
merator will vanish if and only if k is a Weyl aligned null 
direction (WAND) 0, HJ, [H[ , the generalization of 
a PND to higher directions. However, for D > 4, only 
algebraically special spacetimes have WANDs, so gener- 
ically even in vacuum the numerator of /(k) will never 
vanish, and neither will the denominator. 

One may also note that in a coordinate system in which 
the coordinates all have the dimension of length, so that 
the metric components are all dimensionless, the com- 
binations of components of the Riemann tensor and its 
covariant derivatives in both the numerator and denom- 
inator of /(k) have the dimensions of length to the neg- 
ative eighth power, so the ratio /(k), and hence Mi and 
A/2, are dimensionless, pure numerical quantities without 
any units. Generically when they are neither zero nor in- 
finite, their magnitudes give two measures of the strength 
of the gravitational field, increasing for stronger gravita- 
tional fields. However, even in arbitrarily weak gravita- 
tional fields, it is possible to have locations in spacetime 
where the denominator, 25(k), of the fraction /(k), for 
all k, is always anomalously small (or even zero) in com- 
parison with the numerator, [i/(k)] 2 , leading to large (or 
even infinite) values for Mi. As noted above, A/2 has 
greater possibilities of being infinite, since it diverges at 
any location where the denominator vanishes for some k 
for which the numerator does not vanish. Therefore, Mi 
and A2 are not always good measures of the strength of 
the gravitational field, as they can be arbitrarily large or 
even undefined (infinite), but generically they are finite 
and small for weak gravitational fields. 



III. PLANE GRAVITATIONAL WAVES 

Plane gravitational waves are one class of four- 
dimensional vacuum spacetimes which generically give 
nonzero Mi. In fact, Malcolm MacCallum has suggested 
[22| that Type N spacetimes, with all principal null di- 
rections coinciding, may be the only such examples ( "ex- 
ercise for the readers!"). They are also examples of 
VSI spacetimes, with all scalar polynomial curvature in- 
variants vanishing, even though it has long been known 
[12T I22I that they have nonzero Cartan curvature invari- 
ant scalars. 

A general plane gravitational wave that is an exact 
solution of the vacuum Einstein equations has the metric 
[26l [27J that can be put into the form 

ds 2 = -dudv+dx 2 +dy 2 + [(x 2 -y 2 )h'l(u)+2xyti! < (u)]du 2 , 

(3.1) 

where h+(u) and h x (u) are the wave amplitude functions 
for the two polarizations, functions only of the null coor- 
dinate u, and where the prime denotes a derivative with 
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respect to that null coordinate. Geodesies in this metric 
obey the equations 

x" = (xh'l+yti^), 

y" = {-yh'l + xh"). (3.2) 

When the wave amplitudes are small, there exists a 
congruence of geodesies in which, for some period of time, 
the transverse coordinates x and y have the approximate 
form 

x(u) m x a + x a h + (u) + y h x (u), 

y(u) w y - y h + (u) + x a h x (u), (3.3) 

where xq and yo are constants for each geodesic in the 
congruence. These obviously obey Eqs. (|3.2[) when the 
small deviations from (xo,yo) are dropped on the right 
hand sides of Eqs. (|3.2[) . (However, including the devia- 
tions of x and y on the right hand sides leads to long-term 
average accelerations of x and y toward the origin as a 
result of the gravitational attraction of the energy con- 
tained in the gravitational wave.) 

The behavior of this congruence of geodesies illus- 
trates how h+(u) and h x (u) axe the dimensionless am- 
plitudes for the waves. For example, if one has two 
particles in the congruence that are separated purely 
in the re-direction, their proper separation is Aa; ~ 
Axo(l + h + (u)), which oscillates as an oscillating wave 
goes by from the h + polarization. Similarly, two particles 
separated purely in the y direction have proper separa- 
tion Ay ps Ayo(l — h + (u)), which oscillates 180° out 
of phase from the same h + polarization. On the other 
hand, two particles separated along a line bisecting the 
x- and y-axes, say with Axo = Ayo, have proper sepa- 
ration (Ax + Ay)/V% ~ [(Ax + Ay )/V2](l + h x (u)) 
from the other polarization. 

For the exact plane wave metric Eq. (|3 . 1 [) . one can 
readily calculate that the nonzero components of the Ric- 
mann tensor are 

Ruxux Ruyuy h_^_(ll^ 

Ruxuy — ~^~Ruyux — tl x (u*), (3.4) 

plus those related by the symmetries R a bcd = R[ab][cd]- 
Furthermore, the only nonzero covariant derivatives of 
the Riemann tensor are those with respect to u, which are 
the same as the partial derivatives with respect to u. As 
a result, one can quickly calculate that the numerator of 
the ratio /(k) in Eq. ^) is 4(hf + h' x 2 ) 2 {k u f , and the 
denominator is A(h'^' 2 + h' x 2 )(k u f . Therefore, the ratio 
/(k) is independent of the auxiliary null vector k (though 
it is strictly speaking undefined if one takes k = d/dv, the 
propagation null vector of the wave, which has k u = 0, 
but one can easily avoid that special direction for k). 
One immediately gets that both local scalar numerical 
invariants are 

, (hf + h" 2 ) 2 , x 

M=^2= ;,j, 2 + ^,; 2 . (3.5) 



MacCallum has informed me [13] that this value is the 
reciprocal of the squared modulus of a Cartan invariant 
calculated by Jan Aman and put into the CLASSI [28[ 
file wave2.nul on June 9, 1986. So this invariant is not 
new, but its presentation as a simple function defined for 
generic spacetimes may be. 

If both h+(u) and h x (u) are monochromatic sinusoidal 
functions with the same frequency, e.g., 

h+(u) = A + cos (uju + 0+), 

h x (u) = A x cos (luu + <p x ), (3-6) 

then one simply has 

M 1 =M 2 = h 2 + + h 2 x , (3.7) 

the sum of the two w-dependent oscillating squared am- 
plitudes for the two polarizations of the monochromatic 
plane wave. 

For linearized plane waves, the effective average energy 
density of the waves, in a local Lorentz frame in which the 
null coordinate of propagation of the waves is u = t — z, 
is [29j . in units with G = c = 1, the average of (h+ + 
/i'x)/(167r), which for monochromatic waves of frequency 
w is the average of uo 2 (h 2 + + h 2 x )/(16ir) = Af\UJ 2 / (16it) . 
Over one wavelength 27r/w in the z-direction, the aver- 
age energy per transverse area (in the x — y plane) is thus 
the average of Afiu/8 over one wavelength. Then since 
each graviton of the monochromatic wave has energy lo 
in units with K = 1, in Planck units (% = c = G = 1) the 
number of gravitons per wavelength and per transverse 
area is simply the average of A/"i/8. That is, the value of 
M\ = A/*2 averaged over one wavelength of a monochro- 
matic linearized plane gravitational wave is eight times 
the number of gravitons per wavelength and per trans- 
verse area. 

Although the number of gravitons per wavelength and 
per transverse area is a Lorentz-invariant quantity in the 
case of a uniform flux of monochromatic gravitons all 
traveling in the same direction, this number is a nonlocal 
quantity (since the gravitons are not localized). However, 
for classical gravitational fields, the numerical scalar in- 
variants M\ and A/2 defined by Eq. (|2.4p have the advan- 
tage of being purely local. They are well defined for all 
locations where the denominator of f(k) in Eq. (|2.3[) is 
nonzero for all possible auxiliary null vectors k, and also 
when the denominator vanishes for some particular null 
vectors provided that the numerator also vanishes just as 
fast at those null vectors (as in the plane wave example 
above). 

One can have situations in which the numerical scalar 
invariants Afi and A2 are not defined (or are infinite) , if 
the denominator of the ratio /(k) given in Eq. (|2.3[) is 
zero for all possible auxiliary null vectors k (e.g., in dc 
Sitter) for A/i or goes to zero faster than the numerator 
does as one approaches some special null vectors for A/2. 
(The fact that Afi is defined as the minimum of the ratio 
means that if the denominator, which is always nonneg- 
ative, is nonzero for any k, Ai will be finite, though A/2, 
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the maximum of the ratio /(k), would be infinite at any 
location where the denominator is zero while the numera- 
tor remains nonzero for some k, or where both go to zero 
but the denominator goes to zero faster as k approaches 
a value where the denominator vanishes.) 

For the plane waves considered above, this vanishing 
of the denominator for all k is indeed the case if one has 
h'"' 2 + h' x " 2 = 0, where the 4th derivatives of both polar- 
ization amplitudes vanish. If one has a monochromatic 
wave, the numerator would also vanish at the same values 
of u, so that one would get a finite limit of the ratio (ac- 
tually zero in the monochromatic case) if one approached 
it from positions where the denominator did not vanish, 
but if one has a non-monochromatic wave, one can have 
nongeneric situations in which /i"" 2 + h" x " 2 vanishes but 
h"? + h' x remains nonzero, which would make A/i and A/2 
infinite. For example, one could have h + (u) = h x (u) = 
16 A cos (uju) — A cos (2uju) at uju = 2irn for any integer n 
(hypersurfaces of measure zero in the entire spacetime) . 
Or, one could have a very special spacetime in which, say, 
h + {u) = h x (u) = Bu 2 , which would make Afi and A/2 
infinite everywhere. However, for a generic position in a 
generic nonvacuum spacetime of at least four dimensions 
(as well as for at least the plane wave VSI spacetimes), 
A/i and A/2 are well defined and are both nonzero and fi- 
nite. (For generic vacuum spacetime in four dimensions, 
j\fi is zero, and A/2 may be either finite or infinite.) 



IV. SCHWARZSCHILD METRIC 

For the Schwarzschild metric, the numerator, [z/(k)] 2 , 
of /(k) given in Eqs. (|2.1[) and (|2.3[) is fairly straightfor- 
ward to evaluate by hand (though it did take me several 
hours to calculate it, check it, and correct a mistake I 
made originally), but the denominator appeared to be 
beyond my ability to calculate by hand in any reason- 
able time. Therefore, I used GRTensorll Version 1.79 
(R4) [30|| , which also took me several hours, with help 
from James MacKinnon and Andrei Zelnikov, to install it 
and learn the basics of using it, so it would not have been 
worth the effort just for the numerator alone. GRTensor 
enabled me to find that the numerator was exactly as 
calculated by hand (0.07 seconds CPU time) and to cal- 
culate the denominator (44.67 seconds CPU time, over 
600 times as long). 

Initially I did the calculation in the standard orthonor- 
mal frame outside the event horizon, with the spatial part 
of the null vector making an angle a with the radial di- 
rection in a local Lorentz frame of an observer at rest 
with respect to the black hole. After some simplification 
by hand of the result from GRTensor, I got 



/(k) = 



6M 



5r 1 + (1 - 2M) (12 esc 2 a - 15) 



(4.1) 



|32j coordinates, in which the Schwarzschild metric is 
2M 



ds z = - 1 



dv z + 2dvdr + r z d0 z + r z sin 

(4.2) 

Then I used the one-parameter family of null vectors at 
each location, 



1 / 2M A x _ 
k = d v + - 1 x z )dr + -d e 

2 V r I r 



(4.3) 



where x is the one parameter (sufficient for a generic null 
vector for Eq. (|2.3p . because of the spherical symmetry 
of the metric and because of the independence of /(k) 
on the magnitude of the null vector k) . The parameter 
x may take any real value. Outside the horizon, x — 
a/1 — 2M/ r cot (a/2), but the null vector k given by Eq. 
(14. 3|) is well defined for any real x even on and inside the 
horizon, i.e., for all r > 0. 

With this form of the Schwarzschild metric and of the 
null vector, the hand-simplified result of GRTensor for 
the fraction /(k) defined by Eq. (|2.3p is 



/(k) 



_ /2M 

~ V 5r~ 



8 6M\ 



2M 



1 



(4-4) 

One can readily check that outside the horizon, using 
the relation given in the previous paragraph between the 
parameters x and a, that /(k) given by Eq. (|4.ip agrees 
with /(k) given by Eq. (|4~4)) . 

The fraction /(k) given by Eq. (|4.4p . a perfect square 
for the Schwarzschild metric, goes to zero when x is taken 
to infinity (and also when x is taken to zero at all loca- 
tions except at the event horizon, r = 2M). Therefore, 
for all locations in the Schwarzschild metric, 



A/i = min/(k) = 0. 

k 



(4.5) 



This vanishing of A/i appears to be a generic feature of 
four-dimensional vacuum metrics, as we have seen above. 

On the other hand, A/2, the maximum of /(k) over all 
null vectors k, depends on the location. For r > 3M, 
/(k) = 00 at 



4 3M 



± 



'1 M\ (1 5M 
3 ~ ~r~ ) I 3 ~ ~r~ 



so for r > 3M, 



A/2 = max/(k) = 00. 



For r < 3M, the maximum for /(k) is at 

2M 



x 2 = 



1 



(4.6) 



(4.7) 



(4.E 



so for 2M < r < 3M, 



To get /(k) both inside and outside the event horizon, 
it is better to use ingoing Eddington-Finkelstein [29l . l3ll . 



Af 2 = max/(k) = 

k 



3M 



5(3Af - r) 



(4.9) 
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whereas for < r < 2M, 



Af 2 = max/(k) = 

k 



3Af 



1 2 



5(15M - 7r) 



(4.10) 



At the horizon itself (r = 2M), Af 2 = 9/25 = 0.36. 

In summary, for the Schwarzschild metric, A/i , the min- 
imum at each location, over all null vectors k, of /(k) 
given by Eq. (|2.3[) . is zero everywhere, whereas A/2, the 
maximum of /(k), is infinite for r > 3M but is finite 
for r < 3M , given by Eqs. f|4.9|) and (|4.10|) immediately 
above for 2M < r < 3M and < r < 2M respectively. 



V. KERR METRIC 

For the Kerr metric with an orthonormal null basis 
(1, n, m, m) aligned along the eigendirections of the Weyl 
tensor and normalized so that 1 • n = —1 and m • rft = +1 
(and with all other scalar products vanishing), and with 
the one-complex-parameter (z) family of auxiliary null 
vectors at each location being 



k = 1 



zzn + zm + zm, 



(5.1) 



one may calculate by hand and by GRTensor (0.13 sec- 
onds CPU time) that the numerator of /(k) given in Eq. 



(R\ lb R J cld k a k b k c k d ^ 



ajb 



72 2 \H 4 \z\ 



where [32 



^2 



M 



[r — iacost 



(5.2) 



(5.3) 



is the only nonzero Weyl curvature scalar in the Type-D 
Kerr metric. 

My unskilled efforts with GRTensor were unsuccessful 
in calculating the denominator 26(h) of /(k), so that 
I gave up after my computer had run nearly twenty- 
four hours without getting an answer. However, Nicos 
Pelavas [34[ and Malcolm MacCallum [22j kindly calcu- 
lated the denominator for me. Pelavas' output takes up 
about one page of print and is a polynomial in z and z 
with each having powers running from 2 through 6 (25 
terms), and with coefficients functions of M, a, r, and 
9. As I had conjectured, the denominator is proportional 
to \z\ 4 for small \z\, with a coefficient (which I could 
not calculate) of 2 6 3 4 5 2 M~ 4 / 3 \^ 2 \ 10/3 , so for small \z\, 
/ ~ 0.04|M 2 V> 2 | 2/3 M 4 , which has a minimum value of 
A/"i = at z = 0. Thus indeed A/i vanishes for Kerr, as it 
is expected to do for a generic vacuum four-dimensional 
spacetime. 

Malcolm MacCallum [13] used spinor notation to show 
how the denominator, 2<5(k), can be written as the abso- 
lute square of the second covariant derivative of a Weyl 
curvature spinor in a frame with k being one of the null 
basis vectors. He explicitly calculated this complex sec- 
ond derivative and expressed it as a sum of 14 terms, a 



reduction from the 25 terms Pelavas had given me for its 
absolute square. 

From MacCallum's expression, I derived the following 
expression for the fraction /(k): 



/(k) 



2Mp 



5AD 



(5.4) 



where p = \/r 2 + u 2 , u = a cos 9, A = r 2 — 2Mr + a 2 , 
and 



D 




2 P 



3p 2 A 



9 4 , , 2 P 

4 cot a A esc a sin f3 + 4A sin (3 5— 

3 3 /r A 



14A cot a cos j3 



Q 



p 2 A 



(5.5) 



with z = tan (a/2)e 1 ^, where a is the angle between the 
spatial part of k and the radial direction and is the 
angle between the nonradial spatial part of k and the 
9 direction in the orthonormal frame determined by the 
orthonormal null basis, eg = (l+n)/v2, e-r = (1— n)/V2, 
e§ = (m + m)/v2, — (m — m)/(iv2), with the null 
boost freedom in 1 and n and the rotational freedom in m 
and m used to make e? and orthogonal to the Killing 
vectors d/dt and d/d(j>, with 



A = 



a sin 9 



2Mr 



(5.6) 



P = r 4 - 3Mr 3 + 8a 2 r 2 - 6u 2 r 2 + Mu 2 r + 8a 2 u 2 - 7u 4 , 

(5.7) 

and with 



Mu(3r 2 - u 2 ) = Macos9(3r 2 



2 2 
a cos 



9). (5.8) 



The value of A/i , which is zero, is the minimum value of 
/(k), which corresponds to the maximum absolute value 
of D, occurring for \z\ = (a = 0) and for \z\ = 00 
(a = 7r). The value of A/2 is the maximum value of /(k), 
which corresponds to the minimum absolute value of D. 
For r 3> M > a and for 9 not too near either or it (the 
axes of rotation), one can always choose the complex z, 
or equivalently the direction angles a and (3 for the null 
vector k, so that both the real and imaginary parts of 
D are zero. In particular, if cot a cos (3 is fixed to have 
the value Q/(14Ap 2 A) to make the imaginary part of 
D vanish, so long as this value is not too great, one can 
choose j3 (and hence a) to make the real part of D vanish 
as well. 

However, if <5/(14^4p 2 A) is too large, then for all j3 the 
real part of D will be positive when the imaginary part is 
set to zero, so that D can never be zero for any complex 
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z, and /(k) will have a finite maximum value that is A/2. 
One can readily show that asymptotically at large r the 

condition for M 2 to be finite is tan0 ~ N /27/98(M/r), 
so that the region with finite A/2 is a cylinder about the 
axis with asymptotic radius i/27/98Af. 

One can also show that on the axis itself, where A = 0, 
and for r large enough that P > (which in the limit of 
small a is for r > 3M), 

9(r 2 + a 2 ) 3 _ r 2 / f^V _ j^V* 
2 ~ 25a 2 (3r 2 - a 2 ) 2 a 2 V + ^ / V 3rV ' 

(5.9) 

VI. NEAR THE SURFACE OF THE EARTH 

If one approximated the gravitational field near the 
surface of the earth by the Schwarzschild metric, then 
A/i would be zero and A/2 would be infinite. However, 
because of perturbations of the mass distribution from 
spherical symmetry, I would suspect that both the nu- 
merator and the denominator of /(k) given in Eq. (|2.3[) 
would at generic nonvacuum locations be nonzero for all 
directions of the auxiliary null vector k, so that both A/i 
and A/2 would be nonzero and finite. (For vacuum lo- 
cations, after the insights given by 22], I suspect that 
generically A/i would be zero, whereas A/2 could generi- 
cally be either finite or infinite.) 

One might think that a reasonably good estimate for 
/(k) and for its minimum (A/i) and maximum (A/2) val- 
ues at some location would be obtained by taking the 
monopole plus quadrupole contributions to the Newto- 
nian gravitational field of the earth. That might indeed 
be correct for the numerator of /(k) when one is far above 
the surface of the earth and its higher multipole irregu- 
larities, at least if one ignores the contribution from the 
gravitational fields of the sun and moon. However, since 
the denominator has each Riemann tensor with two co- 
variant derivatives, the small-scale variations from local 
mass distributions are likely to give by far the largest 
contribution to the denominator near the surface of the 
earth. They are also likely to prevent the numerator from 
having a minimum (say with the timelike part of the null 
vector fixed) nearly so small as it would be from just the 
quadrupole field of the earth. 

For concreteness, let us fix the magnitude of the aux- 
iliary null vector k so that its timelike part is unity in 
the frame of the surface of the earth (i.e., so that the 
dot product of the null vector with the four-velocity of 
the earth surface is —1). Then in this same frame, the 
Riemann curvature tensor will have typical components 
of the order of magnitude of 

C=U%«1.7xl(T 23 m- 2 , (6.1) 

© c 

where Mq and -R®c 2 are the mass and radius of the earth. 
Nearby objects that distort the spherical gravitational 



field will give distortions in the Riemann tensor of the 
same order of magnitude (more nearly correctly, smaller 
by the ratios of the densities of the objects to the average 
density of the earth, but since these ratios are of the order 
of unity, here I shall ignore them) . 

Therefore, the magnitude of the numerator of /(k) can 
be expected to be of the order of C 4 for generic directions 
of the spacelike part of k. From the results above for the 
Schwarzschild metric, it would vary as sin 8 a, where a 
is the angle between the spatial part of the null vector 
and the radial direction, if the gravitational field were 
precisely spherical and vacuum, but the distortions of 
nearby objects, e.g., in one's office, would presumably 
make the numerator have a magnitude of the order of C 
for generic directions of the spacelike part of k. However, 
in a generic vacuum spacetime, the numerator would be 
expected to vanish as sin a, where a is now the an- 
gle between the spatial part of the null vector and the 
spatial part of the nearest principal null direction. In a 
non- vacuum region, the minimum value of the numerator 
(with the magnitude of k fixed as above) would be ex- 
pected to be of the order of the tracefree part of the Ricci 
tensor, Gp/c 2 ~ Cp/p®, where p is the mass density of 
the location and p® is the average mass density of the 
earth. 

Similarly, the magnitude of the denominator of /(k) 
can be expected to be of the order of C 2 /L 4 for generic 
directions of the spacelike part of k, where L is a typical 
length scale for the smallest nearby object (at a distance 
comparable to its size) that is distorting the Riemann 
tensor of the earth's gravitational field. Therefore, one 
would expect that /(k) would for generic directions of 
the spacelike part of k have magnitudes of the order of 
C 2 L 4 , with its maximum value being of this same order 
(if it is finite) and its minimum value being of the order 
of C 2 L 4 p/p®. If one takes L ~ 1 m, say for the old 
CRT computer monitor some fraction of a meter in front 
of me that I had when I started writing this paper, and 
a location within the air where p/p® ~ 10 -3 , then one 
would expect that, very crudely, 

M ~ KT 26 , A/2 ~ 1CT 23 (6.2) 

(unless A/2 = co) at a typical location in the air, within 
a meter or so of lumpy stuff of roughly the density of 
the earth. These order-of-magnitude estimates exhibit 
the fact that for generic weak nonvacuum gravitational 
fields, I would expect both A/i and A/2 to be very small, 
increasing quadratically with the strength of the field 
(e.g., quadratically with the mass densities and hence 
also with C for a given spatial distribution of the relative 
densities). 

VII. OTHER LOCAL INVARIANTS NONZERO 
FOR VSI SPACETIMES 

The examples of A/i and A/2 defined by Eq. (12. 4|) as 
the minimum and maximum, respectively, of the ratio 
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/(k) given in Eq. (|2.3[) , are of course only two examples 
of an obviously infinite sequence of such examples of lo- 
cal invariants that might be nonzero for VSI spacetimes. 
(Those examples were merely the simplest examples that 
I could readily think of.) All one needs is some contrac- 
tion of Riemann tensors and their covariant derivatives 
with null vectors in both the numerator and denomina- 
tor, so that there are the same number of copies of each 
null vector in both (in order that the ratio not depend 
on the normalization of each null vector, but only on its 
direction, which forms a compact set allowing one to take 
the minimum and the maximum of the ratio as local in- 
variants). One can even have arbitrarily many different 
null vectors in both the numerator and denominator, so 
long as each null vector occurs the same number of times 
in both. Alternatively, one can have different numbers 
of null vectors in the numerator and denominator if one 
normalizes them so that, say, either the numerator or the 
denominator is fixed to be unity. 

Instead of equivalence classes of null vectors related 
by scalar multiplication for the objects to contract with 
the Riemann tensors and their covariant derivatives, one 
could alternatively take equivalence classes of other geo- 
metric objects that are equally contracted with Riemann 
tensor objects in both the numerator and the denomi- 
nator, at least so long as each equivalence class forms a 
compact set so that one can define the maximum and 
minimum of the resulting ratio. 

Of course, it might be difficult to evaluate more com- 
plex examples. Even the simple examples given here ap- 
pear rather intractable much beyond the simple cases 
considered above of a plane gravitational wave and of the 
Schwarzschild metric. However, their conceptual possi- 
bility shows that there are far more local invariants for 
a gravitational field that may be concocted than just the 
polynomial invariants that are normally constructed. 

However, all of these local scalar invariants are just 
special cases of what can be constructed from the Car- 
tan invariants, since they are determined by components 
of the curvature and of its covariant derivatives in special 
bases either partially or totally determined by the curva- 
ture and its derivatives. Therefore, they are not logically 
independent methods of concocting other invariants be- 
yond those that may be obtained from the Cartan invari- 
ants, but some of them might be thought of as somewhat 
simpler concoctions that shortcut the most general form 
of the Cartan procedure. 

It might also be worthwhile being reminded that if one 
does not restrict attention to local invariants, as was done 
here, then there would be an even vastly greater set of 
invariants. For example, even for flat spacetime, if space 
is compactified into a torus, then the periods and angles 
of the torus would be nonlocal invariants. When there is 
curvature, one can readily think of even many more ways 
to construct nonlocal invariants. 



INVARIANTS FOR NULL FLUIDS AND 
ELECTROMAGNETIC FIELDS 



VIII. 



The examples of N\ and A2 defined above by Eq. 
were concocted to be nonzero for generic spacetimes, and 
also for generic plane waves that are VSI spacetimes, 
though apparently A/"i is zero for generic four-dimensional 
vacuum spacetimes. However, one can have VSI space- 
times that are not vacuum but have null fluids, elec- 
tromagnetic waves, dilatons, and/or supergravity fluxes 
[E i, S H i, 0, 1, i , for which in generic circumstances 
one can concoct nonzero local invariants that are simpler 
than A/i and A/jj. 

For example, instead of the ratio /(k) defined by Eq. 
(|2.3p . for a nonvacuum spacetime one could define 



/(k) = 



(8.1) 



in terms of the Ricci tensor R ao (or in terms of its trace- 
free part, S ao = Rab — jjRQab m D dimensions, which 
would give the same result) and then take 



3 = min/(k), A/4 = max/(k). 

k k 



■2) 



Assuming that the covariant derivatives of the Ricci ten- 
sor are suitably nonzero, these invariants may be well 
defined and nonzero even for generic null fluid VSI space- 
times. 

For an electromagnetic field tensor F ao , one might de- 
fine 

= F* a F ib Fi c F jd F k e F kf k a kWkekf 

J\ K )— m 77 1 l.»l.kI.rUl. P [.f 



F\. bc F ld]ef k a k b k c k d k e kf 



and then take 



A/5 = min/(k), A/6 = max/(k). 

k k 



(8.4) 



For simplicity (though more general examples in curved 
spacetime would also be interesting to investigate), let 
us consider the case in which the electromagnetic field 
stress-energy tensor has a negligible effect upon the 
spacetime, which has the nearly flat metric 

ds 2 = -dudv + dx 2 + dy 2 (8.5) 

with u = t — z,v — t + z. Then take the example of a 
plane electromagnetic wave traveling in the z-direction, 
with electromagnetic potential one-form 

A = A x (u)dx + A y (u)dy (8.6) 

and electromagnetic field two-form 

F = dA = A' x (u)du A dx + A' y (u)du A dy. (8.7) 

With u = t — z, this gives E x = B y = —A' x and E y = 
—B x = —Ay. This plane- wave electromagnetic field then 

gives /(k) that is independent of k and has the form 



All/2 1 All/2 



(8.8) 
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If both A x (it) and A y (u) are monochromatic sinusoidal 
functions with the same frequency, e.g., 

A x (u) = a x cos (uju + <j> x ), 

A y (u) — a y cos (wm + (f) y ), (8-9) 

then one has 

M 5 =M 6 = (E 2 x +E 2 y ) 2 /u\ (8.10) 

Since the energy density of the electromagnetic wave 
is (E 2 + £^)/(47r), over one wavelength 2it/uj in the z- 
direction the average energy per transverse area (in the 
x — y plane) is thus the average of V^sW2 over one 
wavelength. Then since each photon of the monochro- 
matic wave has energy u>, the number of photons per 
wavelength and per transverse area is simply the average 
of *JJJ$/2. That is, the value of — y/TJe averaged 
over one wavelength of a monochromatic linearized plane 
electromagnetic wave is twice the number of photons per 
wavelength and per transverse area. 
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IX. CONCLUSIONS 



In this paper I have given simple illustrations of the 
long-known fact that VSI spacetimes (which have all 
scalar invariants vanishing everywhere that are con- 
structed as polynomials in total contractions of the Rie- 
mann tensor and its covariant derivatives) have other 
classes of local scalar invariants that generically do not 
vanish. The examples in this paper are concocted by tak- 
ing the ratio of total contractions of the Riemann tensor 
and its covariant derivatives with the same number of 
copies of an auxiliary null vector in both the numera- 
tor and denominator and then taking the minimum or 
the maximum of the ratio as the direction of the null 
vector is varied over its compact unit sphere on the null 
cone. When these minima or maxima are finite, they are 
well-defined local Lorentz invariants of the gravitational 
field. Examples were given that reduce to the squared 
amplitudes of monochromatic plane gravitational waves 
(and whose average over one wavelength is proportional 
to the number of gravitons per wavelength and per cross- 
sectional area) . The first of these examples seems to van- 
ish for generic vacuum spacetimes in four dimensions. 
Both of these example local scalar invariants were fully 
calculated for the Schwarzschild metric, partially calcu- 
lated for the Kerr metric, and estimated near the surface 
of the earth. Other examples were given for null fluids 
or null electromagnetic fields. 

These examples illustrate the point made by the Car- 
tan constructions [O, EE [3 [H, H3, d, EM H3, Hi| that 
local scalar invariants may be concocted that are nonzero 
for a larger class of spacetimes than one might naively 
think, such as the VSI spacetimes for which all scalar 
polynomial invariants in the curvature and its derivatives 
vanish. 
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